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The suppression of superconductivity by nonmagnetic disorder is investigated systematically 
in the organic superconductor «:-(BEDT-TTF)2Cu(NCS)2. We introduce a nonmagnetic dis- 
order arising from molecule substitution in part with deuterated BEDT-TTF or BMDT-TTF 
for BEDT-TTF molecules and molecular defects introduced by X-ray irradiation. A quanti- 
tative evaluation of the scattering time TdHvA is carried out by de Haas-van Alphen (dH[vA) 
effect measurement. A large reduction in Tc with a linear dependence on 1/rdHvA is found in 
the small-disorder region below 1/rdHvA =i 1 x 10^^ s"^ in both the BMDT-TTF molecule- 
substituted and X-ray- irradiated samples. The observed linear relation between Tc and 1/rdHvA 
is in agreement with the Abrikosov-Gorkov (AG) formula, at least in the small-disorder region. 
This observation is reasonably consistent with the unconventional superconductivity proposed 
thus far for the present organic superconductor. A deviation from the AG formula, however, 
is observed in the large-disorder region above 1/rdHvA — fx 10^^ s~^, which reproduces the 
previous transport study (J. G. Analytis et ai: Phys. Rev. Lett. 96 (2006) 177002). We present 
some interpretations of this deviation from the viewpoints of superconductivity and the inherent 
difficulties in the evaluation of scattering time. 

KEYWORDS: organic superconductor, disorder effect, h:-(BEDT-TTF)2Cu(NCS)2, de Haas-van Alphen 
oscillations, scattering time. X-ray irradiation 



1. Introduction 

Superconductivity in organic charge transfer salts, k- 
(BEDT-TTF)2X, where BEDT-TTF (or ET in short) 
denotes bis(ethylenedithio)tetrathiafulvalene, has been 
investigated extensively. ""^'^ Its remarkable feature is that 
the native quarter-filled band is modified to the effective 
half-filled band by a strong dimer structure consisting of 
two BEDT-TTF molecules. Thus, this family of organic 
conductors has been considered as a typical bandwidth- 
controlled Mott transition system with strongly corre- 
lated electrons.3'4 In the k-(BEDT-TTF)2X system, the 
bandwidth can be controlled by applying physical pres- 
sure and a slight chemical substitution of the molecules, 
which can change the conduction bandwidth W with 
respect to the effective Coulomb repulsion U between 
two electrons on a dimer. A first-order electronic phase 
transition line divides the phase diagram into the su- 
perconducting and antiferromagnetic (AF) Mott insula- 
tor phases.^' The superconductivity near an AF Mott 
insulator has been theoretically expected to possess an 
unconventional d-wave order parameter.^' ^ Experimen- 
tal results on unconventional superconductivity, however, 
have not been elucidated yet.^ 

To clarify the pair-breaking mechanism, the investi- 
gation of the disorder effect on superconductivity is im- 
portant because it can give information on the pairing 
symmetry. ^"^^ There have been several attempts to intro- 
duce disorder into organic superconductors, which arises 
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from molecular defects and impurities by electron, pro- 
ton, and X-ray irradiation^'^^^^ and anion X or donor 
molecule substitution, respectively. In addition, the 
conformational disorder of terminal ethylene groups in 
the BEDT-TTF molecule has been investigated from the 
viewpoint of the cooling speed dependence of electronic 
properties at low temperatures. ^^"^^ In the case of or- 
ganic superconductors, however, there are characteristic 
difficulties in the experimental investigations of the dis- 
order effect. Superconductivity is very sensitive to the 
applied pressure and strain, probably owing to the soft- 
ness of the molecules and molecular lattices. For example, 
the hydrostatic pressure dependence of Tc becomes ap- 
proximately -3 K/kbar in k-(BEDT-TTF)2Cu(NCS)2,^^ 
which is roughly two orders of magnitude larger than that 
of inorganic metals, for example, -0.035 K/kbar for Pb.^'' 
In addition, it has been revealed from Ehrenfest analysis 
for thermal expansion measurements^^ that the uniaxial 
stress/strain effect on Tc is quite large and anisotropic for 
the crystal axes. A large negative uniaxial pressure com- 
ponent of approximately -6.2 K/kbar is found for the in- 
terplane direction, while intraplane components indicate 
positive values of -f3.4 and -f0.14 K/kbar in k-(BEDT- 
TTF)2Cu(NCS)2. Such a large structural effect on super- 
conductivity makes the systematic study of the disorder 
effect difficult in general. 

In this paper, we present a study of the disorder 
effect on the superconductivity of the organic super- 
conductor k-(BEDT-TTF)2Cu(NCS)2, which was syn- 
thesized as the first 10 K class organic supercon- 
ductor.^^ Three different types of disorder are in- 
troduced by the partial substitution of BEDT-TTF 
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molecules with deuterated BEDT-TTF or BMDT- 
TTF molecules, where BMDT-TTF (or MT) denotes 
bis(methylenedithio)tetrathiafulvalene, and molecular 
defects induced by X-ray irradiation into the crystals. 
To discuss the relation between Tc and disorder system- 
atically, we carried out de Haas- van Alphen (dHvA) mea- 
surement for the quantitative evaluation of the disorder 
level using the scattering time, which could be obtained 
more reliably without the complex assumptions needed 
in a transport study. In addition, is determined ther- 
modynamically by measuring the magnetic susceptibil- 
ity. 

From the results, a large reduction in Tc with a lin- 
ear dependence on l/rdHvA, evaluated on the basis of 
the dHvA effect, is found in the small-disorder region 
below 1/rdHvA =i 1 x lO^^ g-i both the BMDT-TTF- 
molecule-substituted and X-ray-irradiated samples. The 
observed linear relation between Tc and l/TdHvA is in 
agreement with the Abrikosov-Gorkov (AG) formula^^ 
at least in the small-disorder region. A possible explana- 
tion for this observation may be the nonmagnetic disor- 
der effect on for the unconventional superconductivity 
discussed so far in this organic superconductor.^^ A de- 
viation from the AG formula, however, is observed in the 
large-disorder region, which reproduces the results of a 
previous transport study. ^'^ We present some interpreta- 
tions of this deviation from the viewpoints of supercon- 
ductivity and the inherent difficulties in the evaluation 
of the scattering time. 

2. Experiment 

Single crystals of k-(/i-BEDT-TTF)2Cu(NCS)2, K-[{h- 
BEDT-TTF)i_^(d-BEDT-TTF)^]2Cu(NCS)2, and k- 
[(/i-BEDT-TTF)i_^(/i-BMDT-TTF):r]2Cu(NCS)2 were 
grown by an electrochemical oxidation method, where h 
and d-BEDT-TTF denote the hydrogenated and deuter- 
ated BEDT-TTF molecules, respectively. Crystals with 
BMDT-TTF molecule substitution could be grown up 
to a substitution ratio oi x — 0.15,^*'^* while crystals 
with partial substitution of the deuterated BEDT-TTF 
molecule were obtained at a: = - 1. The substitution ra- 
tios X of d-BEDT-TTF and /i-BMDT-TTF for /i-BEDT- 
TTF were examined by analyzing the molecular vibra- 
tion of the terminal ethylene group of /i-BEDT-TTF in 
the infrared reflectance spectra. The magnitude of the vi- 
bration modes for the terminal ethylenes should change 
linearly with the number of terminal ethylenes. The 
actual substitution ratio was confirmed by this method 
and the measured ratio x was found to be almost the 
same as the nominal value in the crystal growth. For the 
X-ray irradiation experiments, a crystal of K-(ft,-BEDT- 
TTF)2Cu(NCS)2 was irradiated at 300 K using a nonfil- 
tered tungsten target at 40 kV and 20 mA. The dose rate 
under the present irradiation conditions was expected to 
be about 0.5 MGy/h following a comparison with a previ- 
ous report by Analytis et al}^ The magnetic susceptibil- 
ity measurements using a superconducting quantum in- 
terference device (SQUID) magnetometer at 5 T showed 
no indication of the production of magnetic impurities 
by the molecule substitution and X-ray irradiation. "^^ 

The superconducting transition temperature Tc of the 




r(K) r(K) 

Fig. 1. (Color online) Temperature dependences of (b) the mag- 
netic susceptibility x and (a) the temperature derivative of x in 
K-[(/i-BEDT-TTF)i_:i,(/i-BMDT-TTF):i,]2Cu(NCS)2, and (d) x 
in the X-ray irradiated /t-(/i-BEDT-TTF)2Cu(NCS)2. The mea- 
surements wrere performed at _ff = 3 Oe perpendicular to the 
b-c plane under the zero-field-cooling condition. The definitions 
of Tc and the width ATc are given in the text, (c) Tempera- 
ture dependence of interlayer resistivity of X-ray-irradiated ft- 
(?t-BEDT-TTF)2Cu(NCS)2. Note that the samples in (c) and 
(d) are different samples irradiated under almost the same con- 
ditions. 

samples was determined from the magnetic susceptibility 
at 3 Oe. The definition of Tc is given in the next section 
in comparison to that used in the resistivity measure- 
ments. Magnetic torque measurements were performed 
to detect dHvA oscillations using a precision capaci- 
tance torquemeter.'^^'"^^ The torquemeter with samples 
was directly immersed in liquid '^He or dense '^He gas 
in a refrigerator, which was combined with a 15 T su- 
perconducting magnet at the High Magnetic Field Lab- 
oratory for Superconducting Materials (HFLSM), IMR, 
Tohoku University. In all the measurements of the mag- 
netic susceptibility and dHvA effect, the samples were 
cooled slowly from room temperature to 4.2 K in ap- 
proximately 24 h in order to obtain the same sample 
conditions and exclude any possible disorder effect on the 
conformational order of the terminal ethylene groups of 
BEDT-TTF molecules.i9-2i,24,33 

3. Results and Discussion 

3.1 Disorder effect on superconductivity 

Figures 1(b) and 1(d) respectively show the temper- 
ature dependences of the magnetic susceptibilities 47rx 
of K-[(/i-BEDT-TTF)i_^(/i-BMDT-TTF):r]2Cu(NCS)2 
and the X-ray-irradiated k-(/i-BEDT-TTF)2Cu(NCS)2. 
The demagnetization effect is corrected by using an el- 
lipsoidal approximation for the planar shape of the sam- 
ples. The ambiguity of the shape and size of the samples 
may lead to an error (approximately ±15%) for the full 
superconducting diamagnetic volume — l/47r.'^'* Tc is de- 
termined from the diamagnetic transition curve as the 
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crossing point of tlic interpolation lines from the norm; 
and superconducting regions. 

The Tc of the pristine samples determined here reprc 
duces well the; previous results in which Tc is obtaine 
thermodynamically by specific heat,^^"^* thermal expai 
sion,^^ and magnetization*'^'*^ measurements. In the r< 
sistivity measuremc^nts, in general, the offset temper; 
ture to a zero resistivity seems to be close to the T^. ol 
tained from the thermodynamic measurements, althoug 
the determination of the offset temperature depends o 
the actual measurement accuracy of the sufficiently lo' 
resistivity. Figure 1(c) demonstrates the resistive supei 
conducting transition curves of the X-ray irradiated /■ 
(/i-BEDT-TTF)2Cu(NCS)2. The resistivity is measure 
along the interlayer direction. The changes of the ten 
peratm-e dependences of the resistivity and Tc with th 
irradiation reproduce well the previous results. In con 
parison with the magnetization results in Fig. 1(d), tb 
offset criterion of the resistivity by 0.1% of the normal r( 
sistivity indicates an almost the same Tc determined froi 
the present magnetization measurements and the prev 
ous thermodynamically obtained Tc of a pristine sampli 
Note, therefore, that Tc in the present paper is lower b 
approximately 1 K than that determined as the middl 
point of the resistive transition curve. 

The width of the superconducting transition : 
evaluated to be the half-width of the temper? 
ture derivative curve dx/dT, as shown in Fig. 1(£ 
for K-[(/i-BEDT-TTF)i_^((i-BEDT-TTF)^]2Cu(NCS); 
These conditions and definitions are the same a„ 
those for the measurements of K;-[(/i-BEDT-TTF)i_a;(d- 
BEDT-TTF)^]2Cu(NCS)2 and the X-ray-irradiated k- 
(/i-BEDT-TTF)2Cu(NCS)2. 

The changes in Tc and ATc with the substitution ratio 
X and irradiation time are summarized in Fig. 2. The 
substitution by BMDT-TTF molecules reduces Tc al- 
most linearly with x up to ~ 0.15. The transition width 
increases concurrently with a decrease in Tc. The ob- 
served large suppression of the superconductivity may 
be caused by the disorder effect with the nonmagnetic 
BMDT-TTF molecule substitution. The increase in the 
transition width indicates that the BMDT-TTF molecule 
substitution also causes an inhomogcneity in the super- 
conductivity. On the other hand, the deuterated BEDT- 
TTF molecule substitution does not markedly affect su- 
perconductivity. In the full range of the substitution, the 
largest reduction in Tc by the disorder effect is expected 
to take place at a; = 0.5. The experimental results, how- 
ever, show that Tc increases linearly from .t = to 1 and 
ATc is almost constant. The increase in Tc results from 
a chemical pressure effect on the narrowing of the band- 
width W with respect to U on the dimer.-'^'' Thus, we 
can conclude that the deuterated BEDT-TTF molecule 
substitution does not contribute at all to the suppres- 
sion of the superconductivity. This is also supported by 
the small change in the scattering time with x, which 
is obtained from the dHvA effect discussed in the next 
section. 

The sample irradiated with X-rays shows a decrease in 
Tc with an increase in the irradiation time. The reduction 
in Tc reproduces well the results reported by Analytis 
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Fig. 2. Molecule substitution and X-ray irradiation time 

dependences of Tc and ATc for (a), (b) K-[(?i-BEDT- 
TTF)i_^((i-BEDT-TTF)^]2Cu(NCS)2 and k-[(/i-BEDT- 
TTF)i_:j,(/i-BMDT-TTF):r]2Cu(NCS)2 and (c), (d) 
BEDT-TTF)2Cu(NCS)2 irradiated with X-rays, respectively 
The curves are guides for the eyes. 

et al}^ In addition, the present magnetic susceptibility 
measurements indicate the broadening of ATc in con- 
nection with the reduction in Tc by X-ray irradiation. 
The behaviors of Tc and ATc are similar to those in the 
samples with the BMDT-TTF substitution. This sug- 
gests that X-ray irradiation induces a disorder effect for 
superconductivity owing to the random potential mod- 
ulation caused by defects in molecules. Considering the 
changes in Tc and ATc by the X-ray irradiation in com- 
parison with the BMDT-TTF molecule substitution, 500 
h X-ray irradiation roughly corresponds to the x = 0.05 
- 0.1 substitution of BMDT-TTF molecules. 

To evaluate the disorder effect on quasi-particle scat- 
tering in metals, the residual resistivity po at low tem- 
peratures is a good physical quantity for estimating the 
scattering time Timp by impurity scattering as po = 
m* /ne^Timp, where m* is the effective mass, n is the 
number of the carriers, and e is the elementary elec- 
tric charge. The scattering time Timp obtained from po 
contains mostly the contribution of large-angle scatter- 
ing by impurities. The reduction in the Tc of k-(BEDT- 
TTF)2Cu(NCS)2 irradiated with X-rays has been ex- 
plained by the AG formula with the quasi-particle scat- 
tering time obtained from the residual resistivity.^^ 
There are, however, some disadvantages of using the 
residual resistivity to obtain the scattering time. First, 
the scattering time is obtained indirectly from the re- 
sistivity because it is necessary to determine m* and n. 
Second, it is experimentally difficult to accurately obtain 
the intralayer resistivity for small samples with a large 
anisotropy between intra- and interlayer resistivities. 
In addition, we need to assume the interlayer transfer in- 
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tegral t± in the case that the interlayer resif 
for obtaining Timp.^^'*'* One of the experi] 
ods of quantitatively obtaining the scatterii 
measurement of the dHvA effect. The sc 
can be evaluated directly from dHvA oscil 
dition to the effective mass and the cross-f 
of the Fermi surface. The scattering time 
from the dHvA effect also includes contribui 
angle scattering in addition to that of lar^ 
tering such as elastic impurity scattering. W 
sider different contributions to TdHvA fron 
quantitative consideration of physical prop 
Tc- The dHvA effect of the molecule-substi 
ray-irradiated k-(BEDT-TTF)2Cu(NCS)2 
in the next section. 

3.2 Disorder effect on de Haas-van Alphe 
Figures 3(a) and 3(b) show the magnetic 
and dHvA oscillations of the BMDT-T 
substituted sample with x = 0.03 at 0.44 
of the X-ray-irradiated sample with differe: 
doses, respectively. The magnetic field is ap 
dicular to the conductive h-c plane. At 1 
fields, a large hysteresis in the torque cu 
for the up and down sweeps of the magn 
cause of the magnetic irreversibility owing 
solid states of the type-II superconductor, 
spike structures observed inside the hyste 
suit from the flux jumps appearing in the nc 
vortex state upon sweeping the magnetic 
In the X-ray-irradiatcd sample, as shown 
the torque curves in the superconducting 
a measurement-to-measurement variation 
same sample is measured with increasin 
dose. The reason for the variation may be 
the sample is remounted onto the torquem 
sample is irradiated with X-rays at room temperature. A 
small difference in the setting condition between the sam- 
ple direction and the magnetic field direction nearly per- 
pendicular to the b-c plane sensitively affects the struc- 
ture of the torque curves in the superconducting region. 

The effect of disorder on the magnetic torque is not 
observed in the normal state above approximately 5 T in 
both the BMDT-TTF-molecule-substituted and X-ray- 
irradiated samples. The magnetic torque curves show 
paramagnetic behavior on which the dHvA oscillations 
are superimposed and show no inclusion of magnetic im- 
purities that might be induced by disorder. 

The dHvA oscillations on the torque curves can be 
seen above approximately 5 T even in the vortex liq- 
uid state below the upper critical field i?c2 in the clean 
pristine k-(/i-BEDT-TTF)2Cu(NCS)2.^^ From the dHvA 
oscillations, the cross-sectional area of the Fermi surface, 
the effective mass m*, and the scattering time TdHvA of 
the cyclotron-orbiting electrons on the Fermi surfaces are 
obtained from the oscillation frequency F, and the tem- 
perature and magnetic field dependences of the oscilla- 
tion amplitude, respectively, on the basis of the Lifshitz- 
Kosevich formula. The oscillation part of the magnetic 
torque Alk excluding higher-harmonic components is de- 
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Fig. 3. (Color online) (a) Temperature dependence of mag- 
netic torque curves of K-[(/i-BEDT-TTF)o.97(/i-BMDT- 
TTF)o.o3]2Cu(NCS)2. The magnetic fields are applied 
perpendicular to the b-c plane. The inset indicates the torque 
dHvA oscillations after subtracting the monotonic background 
torque curves, (b) Magnetic torque curves of k-(/i-BEDT- 
TTF)2Cu(NCS)2 irradiated with X-rays with increasing 
irradiation dose. The temperature indicated in parentheses is 
the temperature for each torque measurements. A slightly dif- 
ferent sample setting of the torquemeter and a slight diff'ercnce 
in the magnetic field direction in each measurement after the 
irradiation at room temperature sensitively affect the amplitude 
and structure of torque curves, especially in the superconducting 
region below approximately 5 T. The inset of the figure shows 
the temperature dependence of the magnetic torque curves after 
163 h of X-ray irradiation. 



scribed as 

Alk « H^RtRdRs sin[27r(^ - 7)], (1) 

Id 

for the temperature factor Rt ~ 
[X{m*/ma)T/H]/sinh[X{m*/mQ)T/H], the Dingle 
factor i?D = exp[— A(m*/mo)rD/i?], and the spin 
factor Rs — cos[TTg{mh/mo)/2]. 7 is a phase factor 
that is usually close to 1/2. Here, A = 2Tr^ckB/eh — 
14.69 T/K, mo and are the free electron and band 
masses, respectively, and Td is the Dingle temperature, 
related to the scattering time TdHvA by the relation 
Tu = ft/27rfcBTdHvA- For the torque-dHvA oscillation in 
three and two dimensions, n is 3/2 and 1, respectively. 
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Fig. 4. Magnetic torque dHvA oscillations at T = 0.44 K 
in (upper) pristine k-(/i-BEDT-TTF)2Cu(NCS)2, (left) K,-[{h- 
BEDT-TTF)i_^(rf-BEDT-TTF)^]2Cu(NCS)2, (center) K-[{h- 
BEDT-TTF)i_^(/i-BMDT-TTF)x]2Cu(NCS)2, and (right) k- 
(/i-BEDT-TTF)2Cu(NCS)2 irradiated with X-rays. The mag- 
netic fields are applied perpendicular to the b-c plane. The am- 
plitude of the oscillations is normalized by the value of each 
oscillation at 14 T. 



Here, n = 1 for two dimensions is used. 

Figure 4 shows the dHvA oscihations of pristine 
K-(/i-BEDT-TTF)2Cu(NCS)2 in the upper panel, 
AC- [(/i-BEDT-TTF) i _^ (rf-BEDT-TTF)^] 2Cu(NCS)2 
in the left panel, «;-[(/i-BEDT-TTF)i_a:(/i-BMDT- 
TTF)a;]2Cu(NCS)2 in the center panel, and k-(/i-BEDT- 
TTF)2Cu(NCS)2 irradiated with X-rays in the right 
panel. The results of the oseillation frequency F = 
602 ± 2 T, the effective mass m* = (3.2 ± 0.2)mo, 
and the Dingle temperature Td = 0.48 ± 0.04 K, 
corresponding to the scattering time TdHvA = 2.5 ± 0.2 
ps in the pristine sample, reproduce well the results in 
previous reports. In the present magnetic field 
region below 15 T, we observe oscillations with a single 
frequency that have been attributed to the a-orbit of 
the cylindrical Fermi surface. No magnetic breakdown 
orbit P, which has been observed at higher magnetic 
fields, is detected because of the lower experimental 
magnetic fields than the magnetic breakdown field. 

From the observed dHvA oscillations shown in Fig. 4, 
the oscillation frequency, effective mass, and scattering 
time are obtained on the basis of the Lifshitz-Kosevich 
formula. These parameters are summarized in Fig. 5, 
which are plotted in comparison with the substitution 
content x and irradiation time. 
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Fig. 5. Molecule substitution and X-ray irradiation time de- 
pendences of dHvA frequency F, effective mass ratio 
m*/mo, and scattering time t^hvA obtained by dHvA mea- 
surements in (a), (b), (c) K;-[(/i-BEDT-TTF)i_3,((i-BEDT- 
TTF)x]2Cu(NCS)2 (open circles) and k-[(/i-BEDT-TTF)i_j;(/i- 
BMDT-TTF):,]2Cu(NCS)2 (open squares) and in (d), (e), (f) 
K-(/i-BEDT-TTF)2Cu(NCS)2 irradiated with X-rays (filled cir- 
cles), respectively. The curves are guides for the eyes. 



3.2.1 Substitution effect with deuterated BEDT-TTF 
molecule, K-[(h-BEDT- TTF)i_^ (d-BEDT- 

TTF).,kCu(NCS)2 

The oscillation frequency shows a small monotonic de- 
crease with X from F = 602 ± 2 T for a; = to 597 
± 2 T for X = 1, which is consistent with a previous 
report. This monotonic change in F with x can be ex- 
plained by the chemical pressure effect. The sample with 
a; = 1 is considered to have a narrower bandwidth than 
the sample with x = 0. This is because the full deutera- 
tion of the BEDT-TTF molecule in fi;-(BEDT-TTF)2X 
induces a negative pressure effect of approximately 15-20 
MPa.^^ The change in F with pressure in K-(/i-BEDT- 
TTF)2Cu(NCS)2 has been investigated, and the pressure 
causes F to increase by approximately 0.3%/10 MPa.^^ 
This pressure effect owing to the full deuteration from 
X = to 1 is expected to lead to a decrease of approxi- 
mately 3 - 4 T at approximately F ~ 600 T. The observed 
change in F is consistent with this estimation within the 
accuracy of the measurements of F. 

The effective mass to* of approximately 3.2too does 
not depend on x, as shown in Fig. 5(b), whereas the 
scattering time ranvA shows a concave variation with x 
even though a large sample dependence exists at a; = 
and 0.5. For this deuterated molecule substitution, we 
do not expect a very large disorder effect on the elec- 
tronic properties, as mentioned for superconductivity in 
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the previous section. The randomness induced by the 
substitution, however, may increase the electron scatter- 
ing to some extent. In such a case, the largest disorder 
effect may occur at approximately x = 0.5. Actually, the 
concave feature of TdHvA with x indicates the possibility 
of the phenomenological relation 1/tcxx(1— cc), known 
as the Nordheim relation, in the residual resistivity of 
metal alloys such as copper-gold.^^' The deuterated 
molecule substitution, however, does not cause a large 
spatial modulation by the local impurity potential, which 
results in large-angle scattering Timp for the resistivity 
and superconductivity. Meanwhile, the scattering time 
obtained by the dHvA effect should also be affected by 
the small-angle scattering in addition to the large-angle 
scattering. The observed concave feature of TdHvA with 
X may be explained by the phase shift of the oscillations 
discussed for the weak inhomogeneity in the metal alloy 
system.^^'^^ 

Recently, an alternative model for the decrease in scat- 
tering time obtained from dHvA oscillations has been 
proposed on the basis of the inhomogeneous electronic 
state. ^"'^^ This model shows that the decrease in dHvA 
amplitude is caused by the interference of dHvA oscilla- 
tions with slightly different frequencies, which originate 
from the statistical distribution of the electronic inho- 
mogeneity in the bulk sample. The present observation 
may be an example of such a case because the change in 
F with X is likely to cause microscopic inhomogeneity in 
the electronic states. Neither macroscopic inhomogeneity 
nor a large disorder effect, however, is expected because 
the substitution does not produce any inhomogeneous 
feature in Tc and ATc, as shown in Fig. 2. 

3.2.2 Substitution effect with h-BMDT-TTF 
molecule, K-[(h-BEDT- TTF)i_^ (h-BMDT- 

TTF)^kCu(NCS)2 
The /i-BMDT-TTF molecule substitution has a much 
stronger influence on the dHvA oscillations than the 
rf-BEDT-TTF molecule substitution. The BMDT-TTF 
molecule substitution markedly reduces the oscillation 
amplitude in comparison with the pristine k;-(/i-BEDT- 
TTF)2Cu(NCS)2, and then the magnetic fields in which 
the oscillations become observable shift to higher mag- 
netic fields. 

As is clearly seen in Fig. 5(a), the oscillation frequency 
decreases approximately by 1% with increasing substitu- 
tion ratio X up to 0.15. Some possible causes are consid- 
ered for the decrease in F. First, the negative chemical 
pressure induced by the BMDT-TTF molecule substi- 
tution may be expected to increase the lattice param- 
eters, similar to the case of the d-BEDT-TTF molecule 
substitution. The lattice parameters of the BMDT-TTF- 
molecule-substituted sample, however, show no signifi- 
cant difference from those of the pristine sample. In ad- 
dition, a reduction in unit cell volume by the substitution 
is expected because the BMDT-TTF molecule is smaller 
than the BEDT-TTF molecule. Actually the crystal k- 
(/i-BMDT-TTF)2Cu[N(CN)2]Br has a smaller unit cell 
volume^^ than K-(/i-BEDT-TTF)2Cu[N(CN)2]Br, which 
is the same K-type organic superconductor^^ as the 
present material. In this case, a positive chemical pres- 



sure is expected for the smaller-molecule substitution, 
that is, F should increase with increasing x. This ten- 
dency is in contrast to the observation. 

The other possible origin of the decrease in F is 
the slight modifications of the Brillouin zone and the 
anisotropy of the transfer energies. The slight modifi- 
cations of the Brillouin zone and the anisotropy of the 
transfer energies sensitively lead to a variation in the 
cross-sectional area of the a-orbit because the a-orbit 
is formed in the extended Brillouin zone.''^ In the sister 
compound system mentioned above, the crystal lattice of 
K-(/i-BMDT-TTF)2Cu[N(CN)2]Br is monoclinic,^^ .^^^ile 
that of K-(/i-BEDT-TTF)2Cu[N(CN)2]Br is orthorhom- 
bic.^^ Considering this difference in the crystal lattice, 
the substitution of the smaller BMDT-TTF molecule in 
the present material might not cause a simple isotropic 
contraction of the lattice parameters. 

An alternative idea may be an inhomogeneous modula- 
tion of local electronic states around the substitution site. 
If the substituted BMDT-TTF molecules do not pro- 
vide the normal amount of electron charges, -|-0.5e/donor 
molecule, to anion molecules, the hole charge on the 
BEDT-TTF molecules around the substitution site of the 
BMDT-TTF molecule also varies to maintain the formal 
charge transfer value in the crystal. Such local modula- 
tion of the electronic states may cause a small change 
in the Fermi surface. At the moment, we do not have a 
conclusive explanation for the small change in F with x. 

In contrast to the deuterated BEDT-TTF molecule 
substitution, the scattering time decreases very rapidly 
from TdHvA — 2.5 ps for a; = to 0.5 ps for x = 0.15. 
The decrease in TdHvA is considered to be mostly owing 
to the impurity scattering induced by the substituted 
BMDT-TTF molecules, which may act as point scatter- 
ing centers. With increasing amount of the substitution, 
however, an inhomogeneous distribution of the substi- 
tution may occur in addition to the impurity scattering. 
The inhomogeneity could influence the decrease in TdHvA 
in addition to the impurity scattering. This additional 
damping effect may be the case of the inhomogeneity in 
space discussed above. As we mentioned in the pre- 
vious section for superconductivity, the increase in ATc 
shown in Fig. 2 may support the presence of inhomo- 
geneity. 

Here we estimate the mean free path I because we 
could observe; dHvA oscillations even tliougli a large dis- 
order from impurities and inhomogeneity are induced in 
the crystal. The mean free path I is calculated to be ~ 
25 nm for x — 0.15, where I = fpTdHvA and — 5 x 10^ 
m/s for the a-orbit. I is still larger than the in-plane co- 
herence length ^ ~ 5 - 6 nm.^' This condition of i ^ ^ 
indicates that the superconductivity of the sample with 
the BMDT-TTF substitution should be in a clean limit. 

3.2.3 X-ray irradiation effect in K-(h-BEDT- 
TTF)2Cu(NCS)2 
The X-ray irradiation may cause molecular defects in 
the crystal.^"* The disorder effect of the defects suppresses 
superconductivity, as shown in the previous section. The 
dHvA oscillations are certainly affected by disorder. The 
oscillation frequency F decreases by approximately 1% 
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after 467 h of X-ray irradiation. This change in F may 
originate from the local imbalance of the charge transfer 
between the BEDT-TTF donor and Cu(NCS)2 anion lay- 
ers. It has been expected that the molecular defects will 
be introduced mainly into the anion molecules related to 
the CN ligand for Cu.^^' In the case of the defects in an- 
ion molecules, less charge transfer from BEDT-TTF than 
the normal value of -|-0.5e/molecule may arise locally in 
the crystal owing to a space modulation of the mean 
anionicity of polymeric anions. This local imbalance of 
the charge transfer may induce effective carrier doping^'^ 
and potential modulation in space as disorder^^ in the 
conductive BEDT-TTF layers. It is difficult in princi- 
ple to evaluate the number of defects and the expected 
change in the charge transfer value, which are induced 
by X-ray irradiation. In comparison to the molecule sub- 
stitution with /i-BMDT-TTF, however, disorder caused 
by approximately 500 h of X-ray irradiation seems to 
roughly correspond to that with the substitution con- 
tent X = 0.05 - 0.1 in terms of the change in F. This 
quantitative relation has also been similarly found in the 
behaviors of T^, and ATc as mentioned in the previous 
section. 

The effective mass does not change with the X-ray ir- 
radiation in the same way as the ft,-BMDT-TTF molecule 
substitution. On the other hand, a large decrease in 
TdHvA is observed. This decrease in TdHvA could also be 
explained by the disorder effect caused by the X-ray irra- 
diation, which has been discussed as a local modulation 
of the charge transfer by defects in anion molecules. In 
fact, the reduction in TdHvA is consistent with the changes 
in Tc and F for the /i-BMDT-TTF molecule substitution 
discussed above. 

To understand the effect of the X-ray irradiation on 
the electronic states in detail, further investigation of the 
morphology and nature of molecular defects is necessary 
in the future. 

3.3 Relation between superconductivity and scattering 

time 

In this section, we discuss the relation between super- 
conductivity and the scattering time TdHvA examined in 
the molecule-substituted and X-ray- irradiated k-(BEDT- 
TTF)2Cu(NCS)2. We emphasize here that the disorder 
level in the present crystals is not too large to cause a 
direct influence on superconductivity through the modi- 
fication of the electronic states at the Fermi level. This 
is because the clean limit condition of the superconduc- 
tivity can be applied without significant changes in F 
and m* on the Fermi surface. Note in particular that the 
crystals in which the disorders are introduced are still 
clean for quasi-particles since dHvA oscillations are ob- 
servable. Therefore, we assume in this section that most 
of the decrease in is affected markedly by changing 
the quasi-particle scattering, at least at the low disorder 
level with TdHvA > 1 ps. 

In the case of the deuterated BEDT-TTF molecule 
substitution, almost no disorder effect is observed on 
the superconductivity. The partial substitution acts as 
a parameter for changing the bandwidth, which becomes 
narrower with the substitution. In other words, this re- 



sult microscopically demonstrates that the deuterated 
BEDT-TTF molecule substitution can be regarded as 
a negative chemical pressure applied to the crystal with 
minimal effect on disorder. This fact is also important 
for discussions on the electronic phase separation in k- 
(BEDT-TTF)2Cu[N(CN)2]Br as has been observed near 
the Mott transition, in which the bandwidth is 
controlled by the deuterated molecule substitution. This 
is because the present results ensure that the phase sepa- 
ration in space takes place not by the structural inhomo- 
geneity but by a purely electronic origin in the first-order 
metal-Mott insulator transition. 

Next, we discuss the suppression of superconductivity 
by the /i-BMDT-TTF molecule substitution and X-ray 
irradiation from the viewpoint of nonmagnetic impurity 
scattering. The suppression of superconductivity by dis- 
order in organic superconductors has been investigated 
in several materials using different means of introducing 
disorder. In most cases, however, the amount of dis- 
order has been evaluated from the residual resistivity. 
In the case of the superconductor, one needs to apply a 
magnetic field to suppress the superconductivity or to 
assume the extrapolation function of the temperature 
dependence of the resistivity curve at higher tempera- 
tures above Tc to obtain the residual resistivity at suffi- 
ciently low temperatures. Moreover, from the intralayer 
resistivity in this class of layered organic superconduc- 
tors, it is difficult to obtain a reliable value owing to the 
experimental difficulties in ensuring a homogeneous cur- 
rent distribution. For the interlayer resistivity mentioned 
above, we need to know the interlayer transfer integral 
t_\_, which should be obtained independently from other 
experiments. Moreover, it is usually difficult to determine 
the value accurately. In this situation, one important ad- 
vance made by the present study is that TdHvA can be 
used to evaluate the suppression of superconductivity. To 
the best of the authors' knowledge, this is the first ap- 
proach to evaluating Tc using TdHvA for not only organic 
superconductors but also inorganic superconductors. 

Anderson's theorem indicates that nonmagnetic im- 
purities do not change Tc in s-wave conventional super- 
conductors.^ On the other hand, magnetic impurities re- 
duce Tc because they break the time-reversal symme- 
try of Cooper pairs. The quantitative behavior of Tc 
is given by the AG formula, which was originally used 
to describe the magnetic impurity effect on the reduc- 
tion in Tc in conventional superconductors.^^' ''^ In an 
unconventional case with non s-wave symmetry, however, 
the same formulation has been expected to be applicable 
to the nonmagnetic impurity effect on the reduction in 
Tc.^°~^^'^^ Indeed, a noticeable reduction in Tc by non- 
magnetic impurities has been succcissfuUy applied to the 
AG formula to confirm the unconventional superconduc- 
tivity.^^ It is therefore interesting to examine whether 
the present substitution and X-ray irradiation effects are 
also occur. The transition temperature T^^ affected by 
the impurity scattering is described by the AG formula 
as 
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Fig. 6. Relation between Tc and I/t^hvA for samples of K-{{h- 
BEDT-TTF)i_a,(/i-BMDT-TTF):„]2Cu(NCS)2 (open squares) 
and K-(/i-BEDT-TTF)2Cu(NCS)2 irradiated with X-rays (filled 
circles) . The vertical bars indicate the half- width of the transition 
curves of the magnetic susceptibility as shown in Fig. l(a.). The 
dashed curve is calculated from the AG formula for tag = t,iiivA 
and T^'^ = 10.5 K. The results obtained from the intcrlaycr re- 
sistivity, shown in Fig. 1(c) (open circles), are plotted in relation 
to the interlayer residual resistivity po on the upper axis. The 
vertical bars for the resistivity data indicate the onset and off- 
set temperatures of the resistive transition curves. The relation 
between po and TdHvA is described in the text. 
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where T^q is in an ideally pure system with tag ~^ c« 
and 4'{x) is the digamma function. In the AG formula, 
TAG is supposed to be the large-angle impurity scattering 
in principle. In short, tag — Timp- 

The scattering time for electrons in metals includes 
multiple contributions from different scattering pro- 
cesses.^'''^* Each physical quantity obtained by differ- 
ent experiments has a characteristic contribution from 
different scattering processes to the evaluated scattering 
time. The residual resistivity, for example, is dominated 
mainly by the elastic large-angle scattering by impuri- 
ties. On the other hand, the scattering time obtained 
from the dHvA oscillations includes additional contribu- 
tions from smaller-angle scattering by dislocations, wciak 
potential modulation in space, and so on. Therefore, the 
scattering rate 1 /t^hvA can be described on the basis of 
Matthiessen's rule as 



1/TdHvA = 1/Timp + 1/To 



thcr: 



(3) 



where Tothcr denotes the scattering time excluding the 
impurity scattering. 

Figure 6 shows the relation between the scattering 



time TdHvA and in the samples substituted with h- 
BMDT-TTF (open squares) and irradiated with X-rays 
(filled circles). Tc is determined by magnetic suscepti- 
bility measurements and the vertical bars indicate the 
half- width of the transition curves as shown in Figs. 1(a) 
and 1(b). The dashed curve is the calculated T^*^ based 
on the AG formula of eq. (2) with Tco = 10.5 K, where 
TdHvA is applied for tag in the AG equation. In the same 
figure, we plot the suppression of Tc by the X-ray irradi- 
ation, which is obtained from the interlayer resistivity in 
comparison with the residual resistivity po as shown in 
Fig. 1(c). The residual resistivity at T = is evaluated 
from the quadratic extrapolation of the resistivity curve 
as p{T) = Po + AT"^ below approximately 20 K. The 
vertical bars for the resistivity data show the onset and 
offset temperaturc;s of the resistive transition curves. The 
Tc values evaluated in the previous report by Analytis et 
al}^ correspond to the middle temperature of the bars. 
The correspondence between po (upper transverse axis) 
and 1/TdHvA (lower axis) is considered on the basis of 
the previous arguments on the interlayer resistivity.^^' 
The residual resistivity for interlayer transport is given 

by 



Po 



TTh" 



1 



2e'^m*dt\ 



(4) 



where d is the interlayer spacing and t±_ is the interlayer 
transfer integral. By following eq. (4), the linear relation 
between po and Timp, which corresponds to TdHvA in Fig. 
6, is given by the parameter t±_ = 0.048 meV (0.034 meV) 
while taking d = 1.522 rmi^'' and rn* = 3.2mo (6.5rr?,o). 
Here, m* = 6.5mo in the case for the magnetic break- 
down orbit at higher magnetic fields.^*^'^^'^^'^* Although 
it is a rough estimation, the applied t±_ is in reasonably 
good agreement with the value (0.03 ± 0.01 meV) re- 
ported by Analytis et al}^ and also with an indepen- 
dent determination of t±_ ^ 0.04 meV by an angular- 
dependent magnetotransport experiment. ^^'^'' A small 
difference in the correspondence between po and 1 /TdHvA 
in the X-ray-irradiated samples may be caused by the 
contribution of 1/Tother in eq. (3). This is because po 
tends to be smaller as Po oc (1/TdHvA — 1/Tother) in prac- 
tice. For such quantitative arguments on the transport 
properties, we need further investigation in the future 
while considering the pristine sample-to-sample depen- 
dence. 

The suppression of Tc by both molecule substitution 
and X-ray irradiation is similarly observed for TdHvA- The 
dependence of Tc on 1/ TdHvA is almost linear in the ini- 
tial region below 1/TdHvA — 1 x 10^^ s~^. This behavior 
is consistent with the AG formula of eq. (2). The previ- 
ous results"'^^ for po are also consistent with the present 
results for TdHvA if Timp = TdHvA is assumed. These ob- 
servations indicate that, at least in the small-disorder 
region, the nonmagnetic impurity scattering induced by 
molecule substitution and X-ray irradiation works as the 
pair-breaker expected in the unconventional supercon- 
ductivity. In addition, the good correspondence among 
the scattering time Timp, TdHvA, and tag suggests that 
the contribution of 1/Tother to 1/TdHvA in eq. (3) is rel- 
atively small and thus TdHvA — Timp holds in a weak- 
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disorder region below l/rdHvA — 1 x 10^^ s~^. 

Note, however, that both the present results and the 
previous results by Analytis et al}^ deviate from the AG 
formula in the; larger-disorder region above l/rdHvA — 1 
X 10^^ s~^, in which the decrease in Tc becomes smaller 
with increasing l/rdHvA and po- In addition, such de- 
viation seems to start at a different l/rdHvA in the 
X-ray-irradiated and BMDT-TTF-molecule-substituted 
samples. Looking more closely at the cases of the X-ray 
irradiation, the deviation observed in the present and 
previous studies^^ starts to appear at somewhat differ- 
ent values of l/xdHvA and at equivalent po values. In 
their previous study, Analytis et al}^ proposed a possi- 
ble mixed order parameter with both non-s-wave uncon- 
ventional and s-wave components for the deviation from 
the AG formula. In addition to such an explanation, we 
note other possible causes of the deviation from the AG 
formula as follows. These are based on the experimen- 
tal difficulty in the accurate evaluation of the scatter- 
ing time from either the residual resistivity or quantum 
oscillation measurement. First, for the deviation from 
the AG formula, the insensitive behavior of Tc at larger 
1 / TdHvA may be caused by the increasing contribution of 
1/Tother to 1 /rdHvA as indicated in eq. (3), in compar- 
ison with 1/rimp — 1/tag- The additional contribution 
Tother may result from the effect of the inhomogeneous 
potential modulation in bulk samples discussed above, 
which induces a damping of the; clHvA oscillations but 
does not markedly affect the reduction in T^. Therefore, 
the different effect of Tother on the dHvA oscillations and 
superconductivity may be the origin of the observed weak 
1/TdHvA dependence at larger disorder levels. 

Second, we comment on the relevance of to the eval- 
uation of the scattering time using the AG formula in 
the X-ray irradiation experiments by Analytis et al}^ 
Recently, Sano et al. have reported that the disorder in- 
troduced by the X-ray irradiation of the organic super- 
conductor K-(BEDT-TTF)2Cu[N(CN)2]Br enhanced the 
electron localization, and then the superconductor be- 
came an Anderson-type insulator with increasing irradia- 
tion dose.^* Such a localization effect is enhanced toward 
the Mott critical point of the insulator-metal transition 
with increasing electron correlations. Considering the lo- 
calization effect, the residual resistivity of the strongly 
correlated electron system with disorder does not relate 
simply to the impurity scattering. In the present case, the 
localization effect caused by introducing disorder in n- 
(BEDT-TTF)2Cu(NCS)2 might be weaker than that in 
K-(BEDT-TTF)2Cu[N(CN)2]Br because the former has 
weaker electron correlations owing to its broader band- 
width in comparison with the narrower one in the lat- 
ter.^* Although the weaker effect of the localization is 
expected in the present material, the residual resistivity 
may become larger than the value derived by the impu- 
rity scattering, as described in eq. (4) with increasing 
level of disorder. This additional increase in po niay be 
one of the reasons for the deviation of from the AG 
formula. 

Finally, we mention the reduction in Tc in comparison 
with that in the case of metal superconductors with mag- 
netic transition-metal impurities. A large suppression of 



Tc by a small amount of magnetic impurities is known to 
occur in Al, Zn, and other metal superconductors.^^ For 
example, the rate of reduction of Tc amounts to approx- 
imately —3 K/0.01 at.% Mn in Zn.^® In such a case, the 
rate of reduction can be expected to be proportional to 
J^ScSi/Ep, where J is the exchange coupling between 
the impurity and conduction electron spins, and and 
Si are the spins of the conduction electrons and impuri- 
ties, respectively.^^' ''^ The value of J ~ 1 eV in Zn-Mn is 
approximately two orders of magnitude larger than that 
in organic conductors (J ~ 0.004 — 0.02 eV).^^ Therefore, 
such a magnetic impurity effect on Tc may not account 
for the present large reduction in Tc in organic super- 
conductors even if small impurity spins are induced in 
molecular defects. 

4. Conclusion 

The suppression of superconductivity is investigated in 
the organic superconductor k-(BEDT-TTF)2Cu(NCS)2, 
in which different types of disorder are introduced by 
molecule substitution and X-ray irradiation. To deter- 
mine the relation between superconductivity and dis- 
order, we systematically examined the electronic states 
of «;-(BEDT-TTF)2Cu(NCS)2 with introducing disorder. 
The disorder induces a random electronic potential mod- 
ulation in space with minimal effect on the electronic 
states at the Fermi level. A large reduction in Tc with 
a linear dependence on l/rdHvA evaluated using the 
dHvA effect is found in the small-disorder region be- 
low 1/TdHvA ^ 1 X 10^2 s-i in both the BMDT-TTF 
molecule-substituted and X-ray-irradiated samples. The 
observed linear relation between Tc and I/t^hvA is in 
agreement with the AG formula. This observation is rea- 
sonably consistent with the unconventional superconduc- 
tivity of this organic superconductor that has been sug- 
gested by theoretical and experimental investigations so 
far. 

However, a deviation from the AG formula is observed 
in the large-disorder region above l/rdHvA — 1 x 10^^ 
s^^, which reproduces the results of a previous transport 
study. As reported by Analytis et al.,^^ this deviation 
suggests that superconductivity does not simply obey the 
AG formula on the basis of the single! non-,s-wave order 
parameter, but involves a mixed parameter of s- and non- 
s-waves. Furthermore, in consideration of the reduction 
in Tc with an increase in the scattering time, there are 
experimental difficulties in the accurate evaluation of the 
scattering time, thereby affecting the superconductivity 
inherently. The present dHvA study was motivated by 
the quantitative evaluation of the scattering time, which 
could be obtained more reliably without the complex as- 
sumptions needed in the transport study. Some of the 
difficulties and uncertainties in the evaluation of the scat- 
tering time might be related to the deviation from the 
AG formula. It is clear that the superconductivity of this 
class of organic materials should be investigated further 
to completely understand all the observed phenomena in 
the past 25 years. 
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